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1. Introduction 

Let X(t) denote the position of a diffusing particle at time t, anomalous diffusion is 
characterized by a deviation of the variance (AX 2 (t)) from the linear time dependence 
which follows from the central limit theorem [lj. This deviation is generally associated 
with a change in the Hurst exponent H defined through: 



Normal diffusion with H = 1/2 separates the regime of subdiffusion for < H < 1/2 
from the regime of superdiffusion for H > 1/2. Note that the limiting value H = can 
be associated with a logarithmic subdiffusive behaviour 



(AX\t)) = (X\t))-(X(t)) 2 ~t 



2H 




(AX 2 (t)) ~hr"t. 




f Laboratoire associe au CNRS UMR 7198. 
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Anomalous diffusion is observed in various situational among which one may cite 
turbulent diffusion [3] , transport in amorphous semiconductors [H |5| , diffusion in actin 
networks [6] and living cells 0IE], more generally transport in disordered media as well 
as on fractals [9] [TT]. 

Different models have been proposed to explain anomalous diffusion. Subdiffusive 
behaviour is obtained with the continuous-time random walk [12] [TJ] in which the 
random walker jumps between lattice sites after waiting for some random time r. When 
the waiting time has a broad distribution ("fat tail"), behaving as ip{r) ~ r -i-2i? 
with < H < 1/2, then the variance behaves as in equation (11.11) . with a logarithmic 
correction to the normal behaviour when H = 1/2 [T5]-|17j. 

The problem of anomalous diffusion has been modelled using fractional dynamics, 
a very active field in the last years [2j [18]. In this approach, fractional derivatives 
replace ordinary ones in standard diffusion equations. One may also mention generalized 
Langevin equations which display anomalous diffusion when the noise correlation 
function and the dissipative kernel have long-time tails [T9] . 

Models of random walks with memory may lead to subdiffusive, normal or 
super diffusive behaviour, depending on the value of a memory parameter (201 EI]- 
Superdiffusive behaviour may also result from the presence of long-range step-step 
correlations (see [TO] p 148). 




Exact results have been obtained for diffusion with quenched randomness in 
one dimension (Sinai model) where asymmetric random transition rates lead to the 
logarithmic behaviour of equation (II .2\\ with \x = 4 [22]. When the transition rates follow 
some self-similar aperiodic sequence, the diffusion is normal as long as the fluctuations 
of the environment, characterized by a wandering exponent Q [23], are bounded (Q < 0). 
When the fluctuations are unbounded (Q > 0), logarithmic diffusion is recovered with 
jj, = 2/fi§|. In the marginal case f2 = diffusion is governed by equation (jl.ip . with H 
varying with the amplitude of the aperiodic perturbation [24] . 

The same scenario occurs in the field of critical phenomena for the Hilhorst-van 
Leeuwen (HvL) model [25J. A perturbation decaying as a power of the distance from a 
free surface may be irrelevant, marginal or relevant depending on the value of the decay 
exponent, which plays the same role as fl [25] [33]. It leads to continuously varying 
surface exponents in the marginal case and to a stretched-exponential surface critical 
behaviour when the perturbation is relevant. Alternatively, the perturbation may decay 
in time with similar effects, this has been considered in the case of reaction-diffusion 
processes with a time-dependent reaction rate [311 135] . 

In the present work we introduce a model for anomalous diffusion, where 
diffusion takes place in a space- and time- dependent potential energy landscape. The 
perturbation leads to a rich variety of dynamical behaviours, depending on the values of 
the exponents which govern its space and time dependence. Truly marginal behaviour 
is obtained when the potential is parabolic in space and decays as \jt in time. Then, 

% For a long list of experimental observations of anomalous diffusion see [5] 

§ Note that f2 = 1/2 for a random sequence, which leads to Sinai's result /i = 4as expected. 
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like the surface magnetic and thermal critical exponents of the HvL model [31] [33] , the 
dynamical exponent varies continuously with the perturbation amplitude below some 
threshold value at which the behaviour is logarithmic. 

The paper is organized as follows. The model and the associated Fokker-Planck 
equation are introduced in section 2. Some scaling considerations are given in section 3. 
Exact results about anomalous diffusion in an energy landscape which is parabolic in 
space are presented in section 4. We solve the diffusion equation in section 4.1. We 
study the mean value and variance of X(t) in section 4.2 and the surviving probability 
in section 4.3. This is followed by a discussion in section 5. Some mathematical details 
are relegated to the appendices. 



2. Fokker-Planck equation 

We consider a particle diffusing in one space dimension in a space- and time-dependent 
energy landscape 

X p 

U(X,t) = A — t~ u , -oo<X<+oo, t>0, (2.1) 
V 

with p G Z* and w6K. The associated force is 

dU 



x 



-AX p - x r w (2.2) 
dX v ; 

and the drift velocity is given by the Nernst-Einstein relation 

V x = = -CX'-H-", C=-^4, (2.3) 
k B T k B T v ; 

where D is the diffusion constant, k B is the Boltzmann constant and T the temperature. 
Let P(X, t) X be the probability to find the particle between X and X + X at time t, 
the probability current is the sum of the diffusion and drift contributions 

J x = ~D^ + PV X (2.4) 
and the conservation of probability leads to: 

d_P = JJx =D ^L_ v dP_ p dv x _ 

dt dx dx 2 x ox dx ' y ' } 

Combining equations (I2.3P and (I2.5P one finally obtains the Fokker-Planck equation 
BP f) 2 P BP 

-m =D Bx 2 + c XP ~ l ^ 3x + {p ~ l)CXP " r ^ p ■ (2 - 6) 

3. Scaling considerations 

The diffusion is normal in the absence of a drift, i.e., when C = 0. Then, under a change 
of the length scale by a factor b > 1, one has 

*' = f, t> = L (3.1) 
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and the dynamical exponent z = 1/H = 2 at the unperturbed fixed point. Note that t 
scales like X 2 so that D is dimensionless. 

It follows from a dimensional analysis that Vx/X scales like 1/t. Thus, according 
to Q 

C (X'y- 2 (lf)- u = b z C X v - 2 r" . (3.2) 

Making use of (|3.1|) . one obtains the scaling behaviour of the perturbation amplitude C 
which transforms aqjjl 

C = b z{1 - w)+p - 2 C, z = 2. (3.3) 

When u > p/2 the perturbation is irrelevant since C has a negative scaling dimension 
and decreases under rescaling. Then one expects normal diffusion, at least for small C 
values. 

When uj = p/2 the scaling dimension of C vanishes and the perturbation is marginal. 
The diffusion should then be governed by a fixed line parameterized by C and leading 
to continuously varying exponents if the perturbation is truly marginal. 

When oj < p/2 the perturbation decays sufficiently slowly in time to become 
relevant. C increases under rescaling and the normal diffusion fixed point is unstable. 
The perturbation then drives the system to a new fixed point with a different scaling 
behaviour. 

Let 4>(X,t,C) be some physical quantity with scaling dimension x$. According to 
equations (13 .ip and (I3.3p . (p transforms as 



0' = 0(x\ t', C) = <f (j, 1 b^c\ = b**4>(x, t, 



C) , (3.4) 



at the normal diffusion fixed point. With b — \C\ 2u > one obtains 

0(X, t, C) = X?*4> (Jt-, j-, sgn(C)^) (3.5) 

where 

X c = |C|-Vfr-2") ; tc = X 2 = \C\- 2 ^p- 2 ^ (3.6) 

are characteristic length and time introduced by the perturbation when u ^ p/2. Xc 
and tc diverge for a relevant perturbation and vanish for an irrelevant perturbation 
when C -)■ 0. 

4. Parabolic energy landscape in space 

4-1- Solution of the Fokker-Planck equation 

When p = 2 the Fokker-Planck equation (12.61) becomes 
dP d 2 P dP 

_ = D _ + CXt -._ +Ct -. p . ,4.!, 

|| This scaling behaviour remains valid for non-zero real values of p provided X is replaced by its 
absolute value in (|2.1j) . 
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Let us introduce the new variables 

£ = £(X,t)=Xe^, r = r(t)= /V^t', 

./to 

f(t) = C / (tT" \' = C — 2_ (4.2) 

./to 1 _ w 

and the change of function 

P(X,t)=e^n(£,r), (4.3) 
so that II(£, r) is a probability density for the new variables. One has: 

^=cr" J® n(£, r) + ci r* e 2 ™ 5 + |E , 

^ = e m) ^ = e 3/(t) (4 4) 

<9x ' ax 2 <9£ 2 ' 1 J 



Collecting these expressions in equation ( 14. ip . one obtains the ordinary diffusion 
equation 

dU n 5211 ^ 

w =D w (45) 

For a particle starting from £ = Xo at r = (i.e. X = Xo at t — to for the original 
variables, see (I4.2p ). the solution of equation (14. 5 p is the Gaussian 

Coming back to the original variables leads to 

PPM) - 1 - T f (Jf -^'"H (4.7) 



where 



F(t) = e- 2 ^*V(t) = / exp ( 2C ^- — ) t' . (4. 



to 



1-uj 



This is the solution of equation (14. ip for t > to corresponding to the initial condition 
P(X, to) = 5(X — X ) since f(t) and F(t) both vanish when t — > t . For other initial 
conditions at t = to the solution can be deduced from (14.71) through convolution. 

4-2. Mean value and variance of the position 

The mean value of the position at time t, (X(t)), and the variance (AX 2 (t)) follow from 
equation (14. 71) and, according to (14. 21) . read: 

(X(t)) = X e-M = X exp (-C - - ) , (AX 2 (t)> = 2DF(t) . (4.9) 

Let us now examine how they behave when w is varied. 
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4-2.1. Irrelevant perturbation, u> > 1. The mean position has the following asymptotic 
behaviour when t 3> to : 



(X(t)) ~ X D 



1 + sgn(C) 



(V*) 



w-1 



W- 1 



, Xoo = X exp 



-sgn(C) 



\U)-1 



OJ-l 



(4.10) 



The limiting value is smaller (larger) than Xq and is approached from above (below) 
when C is positive (negative). 

Through the change of variables u = t'/t, F(t) in (I4.8p can be rewritten as: 



F(t) = texp 2C 



t /t 



exp 2Ct 



U 



1-u 



1-U 



u . 



Expanding the exponential and integrating term by term, one obtains 

/ 4.1-u \ 00 1 / \ k rl 



texp 2C 



{Ct 



1 — U)\k 



(t /t)(Ctl 



k=0 
I-.,: - 00 



Jfe!(l-w) 



1 + k{\ - cu) 



1 J 2s k\(l - 

7 fc=0 v 



l + Jfefl -u) 



7 fc=0 v 7 

= t {i + o[(toA) w_1 ] +o(toA)} , 

and .F(t) behaves as: 

' t {1 + 0[(t c /t) w - 1 ]} when 1 < to < 2 
t [1 + O(*o/*)] when u > 2 

Thus, asymptotically we obtain the normal diffusion behaviour 

(AX 2 (t)} ~ 2Dt = 2D X£- , t » t , *c , 

as expected in the presence of an irrelevant perturbation. 



F(t) 



(4.11) 



(4.12) 



(4.13) 



(4.14) 



4-2.2. Marginal perturbation, u> = 1. When u; — )■ 1, /(£) — >■ Cln(t/t ) so that the 
integral giving F(t) in (14. 8[) can be evaluated exactly: 



t 



5 mF(t) = t-/JO^f= 2C+1 



i-(toA) 



2C+1 



Making use of these results in (14. 9p . one obtains: 

(x(t)) = x (t/t y c , 



(AX 2 (t)> 



2Dt 
2C+1 



1 - (t /t) 



2C+1 



when C > -1/2 
2Z>tln(t/t ) when C = -1/2 

(t/t ) 2|C| [l - (V*) 2|Chl l when C < -1/2 



L 2|C| 



(4.15) 



(4.16) 



The mean position is algebraic, with a continuously varying exponent for any non-zero 
value of C. It decreases to zero when the potential is attractive (C > 0) and grows to 
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Figure 1. Anomalous diffusion from the initial position Xq = 100 at to = 1 in a 
marginal attractive parabolic potential decaying as The probability density (top) 
and the reduced potential energy (bottom) are shown for C = 1/2 and D = 1/2 at 
times t = 100 (blue), 200 (green), 400 (red). The variance has a normal linear time 
dependence whereas the mean position decays as 1/ y/i 



0.02 




Figure 2. Anomalous diffusion from the initial position Xq = 4 at to — 1 in a marginal 
repulsive parabolic potential decaying as t . The probability density (top) and the 
reduced potential energy (bottom) are shown for the threshold value C — —1/2 and 
D = 1/2 at time t = 100 (blue), 200 (green), 400 (red). Here the variance has a 
logarithmic correction and grows as t In t whereas the mean position increases as . 



infinity when the potential is repulsive (C < 0). When C > —1/2, the leading behaviour 
of the variance is normal with a dynamical exponent z = 2 ((AX 2 (t)) ~ t 2 / 2 ) whereas 
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Figure 3. Anomalous diffusion from the initial position Xq = 1 at to = 1 in a marginal 
repulsive parabolic potential decaying as i -1 . The probability density (top) and the 
reduced potential energy (bottom) are shown for C = — 1 and D = 1/2 at times t = 100 
(blue), 200 (green), 400 (red). For this value of C < — 1/2 the mean position increases 
as t and the variance scales like (X(t)) 2 oc t 2 . 



the amplitude is C-dependent. The marginal behaviour shows up in the correction-to- 
scaling. A logarithmic correction to the normal diffusion appears at C = —1/2. Below 
this threshold value the dynamical exponent is continuously varying, as expected in the 
presence of a truly marginal perturbation. Since z{C) = 1/\C\ < 2, a superdiffusive 
behaviour is obtained when C < —1/2. 

The time evolution of the marginal probability density and the corresponding 
reduced potential energy are shown for C = 1/2 in figure 1, at the threshold value 
C = —1/2 in figure 2 and below the threshold, for C = —1, in figure 3. Note the change 
in the length scales of the three figures. 



4-2.3. Relevant perturbation, < u < 1. According to (14. 9 j) the mean position behaves 
as: 



(X(t)) =X exp 



-sgn(C) 



if/tcf- 



(to/to) 1 - 



1-u 



(4.17) 



Thus the mean position vanishes (diverg 6s) clS db stretched-exponential function of the 
time when the potential is attractive (repulsive). 

When C > the main contribution to F(t) in (14. 8p comes from the vicinity of 
t' = t. Expanding f(t') near t to second order in u = t' — t, one obtains 



F(t) 



o 

to-t 
X 



I _ e -2Ct- u (t-t ) 



2w-l 



' ■ 2Ct-" 8C 2 

| 2 - e -2Ct-<"(t-*0 [ 4 ^ 2 t- 2 -(t - t ) 2 + ACt~ w {t - t ) + 2] } + • ■ 
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2C 



{i+o[(w*) 1_w ]} 



(4.18) 



when t 3> to and tc, in agreement with the exact results of appendix A in the same 
limit. Thus the variance is given by: 



(AX 2 (t)) ~ ~ DX^ 



t > to, t c 



When C < we can split the integral giving r(t) in two parts such that 
r(t) = r(oo) - 5r(t) , r M = lim r(t) = Jfo = t ) , 5r(t) = /(^ = t) 
with, according to (14. 2ft : 



t— >oo 



/(tl 



cxp 



-2IC 



1 -u 



t' 



(4.19) 



(4.20) 



(4.21) 



This integral is studied in appendix B. 
When t > t c , 6r(t) is given by (IB. 5ft : 

t 



5r(t) 



2\C\ 



exp -2\C\ 



J.1 — £J J.1 — CJ 

^ ^0 

1-w 



(l + O^/t) 1 -]} 



When t > t tc, Too is given by ( IB. 51) . too, and reads: 



if 



2 CI 



{1 + O [(tc/to) 1 -"] } 



Thus, using (14. 8ft . we obtain: 

t 



Fit) 



2C 



exp 2|C| 



J.1 — CJ J.1 — W 

1- w 



2C 



(4.22) 



(4.23) 



(4.24) 



This expression is valid for large values of |C| such that tc <C t and tc <C to- It agrees 
with the exact results in appendix A in the same limit. 
When t > tc > t , ^ is given by (jBTfll 



[1 + O(t /tc)] 



7"nn ft/. 



2C 



1-w 



r 



1 -UJ 



(4.25) 



so that: 



^(0 ~ a uj ^ ex P 



2|C|t 



1-w 



l-UJ 



2ICI 



(4.26) 



once more in agreement with the results of appendix A. 

Collecting these results in (14. 9p gives the long-time behaviour of the variance 

D kexp(2|C| — 



(AX 2 (t)> 



|c| L 



J.1 — UJ j.1 — Ul 



id 



t>t >t c , (4.27) 



and 



{Ax2{t)) ~W\ 



a u t c exp 



2|C|t 



l-O) 



1-W 



, t > tc > t . 



(4.28) 



In terms of scaled variables, we obtain 



<AX 2 (t)> 



D 



cxp 



^"-(to/tc) 1 -" 

l-UJ 



t 

T c 



t>t >f c , (4.29) 
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Table 1. Variation with ui and C of the long-time behaviour (t 3> tc ^> to) of the 
mean position (X(t)), the variance (Al 2 (t)) and the surviving probability S(t) for 
the diffusion in a parabolic potential decaying as t~ w . 

w C <X(i)) (AI 2 (t)) 5(i) 

>1 Any X cxp(-^|>-r) 2Dt cx i" 1 / 2 

1 >-l/2 X (tAo)- C Mi 1 cxi^ 1 / 2 

1 -1/2 X ox /t/h 2Dt\n(t/t ) cx (lnt)- 1 / 2 

1 <-l/2 Xo(tAo) |C| ^pT*o(i/io) 2|C| 5 00 +a<- 2 l c l+ 1 

<1 >0 Ar cxp(-C^) %t u ocr^e 

<1 <0 X exp(|qi^) a u ^t c exp(2|C|^) S^+^e 



2kV 



and 



(AX 2 (t)> 



f2(t/tc, 
a^exp — 

V 1 — u 



l-O) 



t>t c >to- (4.30) 



With an attractive potential (C > 0), the variance is subdiffusive with an in- 
dependent exponent whereas it is super diffusive with a stretched-exponential behaviour 
when the potential is repulsive (C < 0). 



4-3. Surviving probability 

In order to study the surviving probability in the half-space X > of a particle starting 
from Xq > at to we introduce an absorbing boundary condition at X = 0, P(0, t) = 0. 
The solution of the Fokker-Planck equation ( 14. ip with a parabolic potential is obtained 
using the image method and reads 



P(X,t) 



^AuDF(t) 



cxp 



(X-X e 



cxp 



(X+X e-/W) 2 



(4.31) 



4DF(*) / V ADF(t) 

according to ( 14. 7p . Let J-{t) t give the probability that the diffusing particle is absorbed 
at X = between t and t + t, the first-passage probability density is given by [36] 



J r (*) = -Jx(0,t) 



1J tttz X=0 



X e 2 f® 



: exp 



•f„ 2 



(4.32) 



dX '"- u ^ttDt 3 ^) V 4Pr(t); 
Note that the drift term does not contribute directly to Jx(0,t) since Vx(0,i) vanishes 



when p = 2 (see equation ( 12. 3D ). The surviving probability can be written as: 

/*oo l*t 

S[t)= / P(X,t)X=l- / F{t')\' . 

In the last integral the change of variable u = Xq / yjADr{t') leads to 

2 



J-(t') t' 

to V 71 " Jx /y/4Dr(t) 



e" n u, 



(4.33) 



(4.34) 
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so that 



2 rXo/y/4Dr(t) ( y \ 



where erf(x) is the error function (see [37] p 297) with the series expansion: 

o 00 2n+l 

erf<«) = ^(-1)"^^. (4.36) 

v n=0 v y 

In the following we apply these results in the different regimes. 

4-3.1. Irrelevant perturbation, lo > 1. Inserting (14.131) into (14.81) . to leading order, one 
obtains 

r(t) ~texp( 2Ct \ } , (4.37) 



u-1 



and (I4.35P gives: 



S(t) ~ erf 



X ( Ct 
exp' 



1-0) \ n 



(4.38) 



_\/4Dt V w-1 

As expected for an irrelevant perturbation, the time dependence is the same as in the 
solution of the unperturbed problem (36] which coincides with ( I4.38P when (7 = 0. 
In the limit t — > oo, the surviving probability Soo vanishes, even in the presence of a 
repulsive potential. A t -1 / 2 long-time behaviour is given by the first term in the series 
expansion (14.361) : 

S(t) -4^^\/t expf-sgn^M^^j t ^ (439) 

7rL> A c y i V u — 1 J 

4-3.2. Marginal perturbation ui — 1. From equations (14. 8p and f 14 . 1 5 j) we deduce: 

{ w^r(t/t ) 2C+1 when C> -1/2 

T(t)=t ° 2C+1 "1 toln(t/to) when C = -1/2 . (4.40) 

( 2jfe[l-(V*) 2|Chl ] when C< -1/2 

The surviving probability is given by: 

erf (v*VwJP^) whenC^ -1/2 
S(t)=< / J \ • (4-41) 



erf . x ° when C = -1/2 

\y/4Dt \n{t/t )' 



In the long-time limit, t 3> to, one obtains: 

fX 0V /{2C + l)/(7rDt ) (t A) C+1/2 when C> -1/2 

X /^nDt \n(t/t ) when C = -1/2 . ( 4 . 42 ) 

erf (X 0V /(2|C| -1)7(4^0) [l + |(t A) 2|Chl ] ) when C < -1/2 



Anomalous diffusion 



12 




Figure 4. Time dependence of the surviving probability of a particle starting from 
Xq = 1 at to = 1 in the marginal case (lu = 1) with C = 1/2 (blue), C = (green), 
C = —1/2 (orange) and C = — 1 (red). S decays very slowly as (lni) -1 / 2 when 
C = —1/2. The decay is algebraic when C ^ —1/2, with a non- vanishing asymptotic 
value when C < —1/2. 



Soo vanishes when C > —1/2. Only when C < —1/2, i.e., when the repulsive force is 
strong enough, does the surviving probability have a non-vanishing limiting value: 



5 



erf X n 



'2\C\ - 1 



C < -1/2. 



(4.43) 



4£>t 

A first-order expansion of the error function gives the deviation from this limiting value: 

2\C\-1 



6S(t) = S(t) -S^-Xq 



2\C\ 



4nDt 



exp 



(2|C|-1)X, 



ADtr 



C<-l/2. (4.44) 



Thus 5S(t) exhibits an algebraic behaviour with different exponents above (t^ ^ 1 ^ 2 ) and 
below (t~ 2 ' c '' +1 ) the critical value C = —1/2 where the decay is logarithmic [(lni) -1 ' 2 )]. 

The time dependence of the surviving probability in the marginal case for different 
values of C is shown in figure 4. 



4-3.3. Relevant perturbation, < cu < 1. When C > and t 3> t and tc, 
equations ( 14. 8 j) and ( 14.18P lead to: 



t w (2Ct x - u 

r(t) "2C exp iT^T 

The surviving probability behaves as: 



(4.45) 
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Figure 5. Time dependence of the surviving probability of a particle starting from 
Xa = 1 at to = 1 in the relevant case (to = 1/2) with C = 1 (blue), C = (green) and 
C = —1/2 (red). When C — 0, the diffusion is normal and the surviving probability 
decays asymptotically as £ -1 / 2 . There is a stretched-exponential decay when C / 0, 
with a non-vanishing asymptotic value when C < 0. 



The asymptotic decay to Soo = has the following stretched-exponential form: 



2 Xq (tc 
T 



2 / {t/tcf^ 



cxp 



nDXc \t J r \ 1 - w 
When C < and t > £ > t c equations fl4~22l and fl423|) give: 



"(*) 



if 



t u 



exp 



2\C\t 



1-w 



2|C] 2|C| A V 
Hence, the surviving probability behaves as: 



~ erf X 



\C\_ 

2Dt U r 



1 ft 



2 Vt 



exp 



2|C|t 



1-0) 



1-W 



It follows that: 



erf I 



1 Ac 



o>/2^ 



(4.47) 



(4.48) 



(4.49) 



(4.50) 



A first-order expansion of the error function gives the deviation from 5^, which decays 
as a stretched exponential: 



8S{t) 



1 X ( t c 



V2nD X c V to 



w/2 

exp 



2Dlx c 



to 



iVex f ggM^! 
v exp v i-w 



When C < and t > tc > to equations fl4~22]) and (14T25]) give: 



( /-i 



2lClt 



(4.51) 
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Comparing with equation 04.48 p . we note that the only change is the substitution of 
a^tc for £q. Hence, the asymptotic surviving probability now reads 

Soo ~ erf ( - 1 . (4.53) 

and the deviation from iSoo displays the same stretched-exponential decay as above: 



5S(t) ~ . = ^-exp 



1 /^o x2 



2a aJ L> V X, 



t V / 2{t/t 



The time dependence of the surviving probability in the presence of a relevant 
perturbation (u = 1/2) is shown in figure 5 for different values of C. 

5. Discussion 

In this work we have shown that diffusion may become anomalous under the influence 
a potential U(X, t) varying as X p in space and decaying as t~ u in time. Scaling 
considerations indicates that the perturbation is relevant (irrelevant) when u < p/2 
(oj > p/2). Then the perturbation introduces a characteristic length Xc and a 
characteristic time tc = Xq, varying as powers of the perturbation amplitude C. The 
scaling functions can be expressed in terms of the scaled variables X/Xc and t/tc- 
When bj = p/2, C is dimensionless and the perturbation is marginal. 

Actually, truly marginal behaviour with C-dependent exponents is only obtained 
for p = 2. A great variety of scaling behaviours then shows up when u and C are varied. 
Our results are collected in table 1. 

The peculiarity of the parabolic potential appears clearly in equation (13.31) . For 
p = 2, the perturbation amplitude C scales as 

C = jf(i-*>)C, (5.1) 

and thus is marginal when oo — 1. Since the dynamical exponent enters as a factor of 
(1 — oS) into the scaling dimension of C, the perturbation remains marginal even when 
z varies with C . This is no longer true for any other value of p where a variation of 
z would immediately lead to a non-vanishing scaling dimension of the perturbation^. 
This may be verified with a linear potential {jp — 1 in equation (12. 61) ) where the problem 
can be solved with the change of variables 

i = X + C— °—, r = t-t , (5.2) 

1 — UJ 



leading to 



(X(t)) =X -C l ° , (AX 2 (t)) = 2D(t - 1 ) , (5.3) 

1 — CO 



% In this respect the HvL model is different, since there the marginal behaviour affects the surface 
magnetic exponents which are continuously varying, whereas it is governed by invariant bulk 
exponents [261127] , 
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for t > t . Here nothing special occurs at oj = 1/2 when the scaling dimension of the 
perturbation at the unperturbed fixed point vanishes. 

With the marginal parabolic potential one may notice that even if the variance 
behaves normally when C > —1/2 there are traces of the marginal behaviour showing 
up in the correction-to-scaling of the variance as well as in the mean position which grows 
or decays as t~ c , depending on the sign of C. The marginal behaviour is apparent, for 
the variance and for the mean position as well, when C < —1/2. Actually (AX 2 (t)) 
scales like (X(t)) 2 when C < -1/2. 

In the relevant situation, u < 1, the mean position always displays a stretched- 
exponential behaviour. The variance behaves quite differently for attractive and 
repulsive potentials as in the HvL model [30]. In the first case, C > 0, (AX 2 (t)) 
increases as t u whereas the growth is of the stretched-exponential type with the same 
scaling as (X(t)) 2 when C < 0. 

The time evolution of the surviving probability with (i) in the marginal situation, 
a logarithmic behaviour at C — —1/2, an algebraic decay to Soo with C-dependent 
exponents and a non- vanishing value of when C < —1/2 and (ii) in the relevant 
situation, a stretched-exponential decay to iSoo and, when C < 0, a non- vanishing value 
of Soo, is the same as the temperature evolution of the surface magnetization m s (T c — T) 
of the HvL model in the ordered phase [31j [32] with 1/t playing the role of T c — T. 

Note that, for the transverse-field Ising model in one dimension, similar 
correspondences have been noticed before between the finite-size properties of the 
surviving probability of a random walk in a random or aperiodic environment and 
the finite-size behaviour of the surface magnetization of the Ising model in the same 
environment [38j-[40j. 

Finally let us mention that the diffusion problem can be solved in the same way in 
higher dimensions for a time- dependent spherically symmetric potential. 

Appendix A. Calculation of F(t) when u = and 1/2 

When to = the function F(t) defined in equation (14.81) yields: 




(A.l) 



When uj 



1/2 we have: 




to 



(A.2) 



Thus, using 




(A.3) 



one obtains: 




(A.4) 
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The approximate results in equations (14.181) . (14.241) and (I4.26P are all in accordance with 
these exact results when the appropriate limits are taken. 



Appendix B. Calculation of I(t\ 



According to (14.21 j) we have: 



1(h) = exp 



1-u 



exp 



2\C\ (t 



/\1— O) 



1-u 



t'. 



(B.l) 



With the change of variable u(t') = 2|C|(t') 1 w / (1 — u), I(ti) can be rewritten as 



where 



i(h 



T(a, x) 



2\C\ 



o) \ T 



exp 



2\C\t 



1-0) 



1-U 



r 



1 2\C\t\ 



1-0) 



1 — U 1 — U 



u a - l e~ u n 



(B.2) 



(B.3) 



is the incomplete gamma function. When t 1 ^> t c we have u(tf) ^> 1. Thus, in this 
limit, the asymptotic expansion (see [37] p 263) 



V(a,x) ~x a -\- x 



a-1 (a-l)(a-2) 
1 + + - ^ '- + 



leads to: 



/(fi)~^exp(-2|C| 



2\C 



x 



j.1— 0) j.1— U) 



1-U 



1 + 



u 



2\C\t[ 



1-0) 



+ 



(x -> oo) (B.4) 



(B.5) 



When ti t c then u(t%) 1 and the series expansion (see [37J p 261-262) 

oo 

T(a, x) = T(a) - x a ^ ,. v ,7,„, (B.6) 



n=0 



(a + n)n! 



can be used to write: 

f(l-uY\~ 



I(tl 



2\C\ 



exp 



2|cn 



1-0) 





1-w 



1 — u I \{1 — u 



2\C\ 



(B.7) 
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